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Abstract 



Recent experiments on Sr2Ru04 have indicated the presence of hnear 
nodes in the superconducting order parameter. Among the possible spin 
triplet states, 2D p-wave superconductivity with symmetry appears to 
be inconsistent with the experiments, whereas the 2D f-wave order parameter 
with Big X Eu symmetry turns out to be a more likely candidate. Here it 
is shown how the quasi-particle bound state wave functions around a single 
impurity provide a clear signature of the symmetry of the underlying super- 
conducting order parameter. 
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I. INTRODUCTION 



In recent years, we have experienced a rapid improvement in the quahty of Sr2Ru04 single 
crystals. This development has lead to a renewed discussion regarding the symmetry of the 
underlying superconductiving order parameter. The data from specific heat measurements 
|T|], NMR results on Ti B , and magnetic penetration depth experiments B in very pure 
single crystals with Tc ^ 1.5K suggest a linear nodal structure in the superconducting order 
parameter, similar to the dx2_y2-wave superconductors Hence, the initially proposed 
fully gapped 2D p-wave superconductivity with symmetry [^] appears to be inconsistent 
with the experimentally observed nodal structure. On the other hand, the spin triplet 
nature of the superconductivity in Sr2Ru04 has clearly been established by muon rotation 
experiments @] which probe a spontaneous spin polarization, and by the fiat Knight shift 
seen in NMR 0. 

Motivated by these recent experiments, several spin triplet order parameters with linear 
nodes have been proposed. [§] In particular, from a study of the thermal conductivity tensor 
in a planar magnetic field it has been suggested that the 2D f-wave order parameter with 
Big X Eu symmetry is a likely candidate. [^] Furthermore, recent measurements of the angular 
dependence of the upper critical field have detected a in-plane anisotropy consistent with a 
Big X Eu component to the order parameter. |Tn 

The objective of this paper is to calculate the quasi-particle bound state wave function 
around non-magnetic impurities in unconventional superconductors with the anisotropic 
order parameters which have been proposed for Sr2Ru04. It has recently been observed in 
a series of STM imaging experiments that a Zn-impurity in the high-Tc compound Bi2212 
gives rise to such a bound state, leading to a distinct four-fold symmetric pattern in the 
local tunneling conductance around the impurity. From the theoretical side it was shown 
that this type of bound state wave function can be interpreted in terms of the solutions of 



the Bogoliubov-de Gennes equations for dx2_y2-wave superconductors. [|T2|,|T3[. Here we will 
perform an analogous analysis for the proposed 2D p-wave and f-wave order parameters. 
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II. BOUND STATE WAVE FUNCTIONS 



A. 2D p-wave superconductors 

For the initially proposed fully gapped odd-parity 2D p-wave superconductivity with an 
order parameter A(k) = A^exp (±i0) the Bogoliubov - de Gennes equations have already 
been worked out. [ P3| , P^ In particular, when A(r) = A = const. [|T1|, one obtains 



Eu{r)= (^-^-^-V{r)ju{r)+p-p'A{id,-dy)v{r), (1) 

Ev{r) = - (^-^ -11- V{v)^ v{v) + pp^/\{id^ + dy)u{v), (2) 

where /i is the chemical potential, and the impurity potential, centered at the site r = 0, is 
approximated by V^(r) = a5^(r). For simplicity, only the 2D system is considered here. 
A simple variational solution of Eqs. (1) and (2) is found to be 

u{r) = Aexp i--fr)Jo{pFr), (3) 
v{r) = Aaexp {—'yr)Ji{ppr) exp (4) 

where pp is the Fermi momentum, Ji{z) are Bessel functions of the first kind, and A is a 
global normalization factor. By inserting these variational wave functions u{r) and v{r) into 



the Bogohubov - de Gennes equations, it follows that E = —(V/2) ± y{K — Vy + A^ and 
a = A/{E + K) ~ 1, where the kinetic K and the potential V contributions to the energy 
are defined by 

drr (9,. exp {-'jr)Ji{pFr)f + (/exp {-^r)Ji{pFr) /rf 



K 



2m drr {exp {—'yr)Ji{pFr) f ' mpp 

The squares of the wave functions u(r) and v{y) can be observed by scanning tunnel- 



^3 

/X ^ (5) 



ing microscopy. |[TT| The tunneling current, /(r, F) oc / dEAs{r,E)AN{r,E + eF), is a 
convolution of the local one-particle spectral function of the normal-state tip, An{y^E) = 
J2k^{E — Ek), and that of the superconducting sample, As{r,E) = J2k[W{''^)\'^^{E — E^) + 



\v{r)\'^S{E + Ek)]. The differential tunneling conductance is thus obtained by taking the 
partial derivative of J(r, V) with respect to the applied voltage V, 



^^,r, V) (X sech^ (^^j I^WI^ + sech^ l^^j '^^1 ■ (7) 

At small temperatures, the local tunneling conductance around the impurity site is domi- 
nated by |'u(r)p for a fixed binding energy Eq and by |f (r)p for —Eq. 

In Fig. 1(a) and (b), we show |M(r)p and |f(r)p respectively. The Fermi wave vector 
Pf — 2.7 /a was chosen to be consistent with band structure calculations. Jl^ The patterns 



described by the squares of these wave functions are concentric circles without a trace of 
four- fold symmetry. If the underlying superconducting order parameter is indeed A(k) = 
Az exp (±i0) the patterns of the local tunneling current around non-magnetic impurities 
should thus be featureless along the azimuthal direction in the 2D plane. 

B. 2D f-wave superconductors 

Let us now consider 2D f-wave superconductors with an order parameter A(k) = 
cos (20) exp (±i0). Following the above procedure, the corresponding Bogoliubov - de 
Gennes equations are then given by 

Eu{r) = f-^ - - Vir)] u{r) + p-/ A{d', - d'^){zd, - a,)i;(r), (8) 
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Ev{r) = - (^-^ - /i - V{r)j v{r) + pfA{dl - aj)(z9. + dy)u{v). (9) 

Bound state solutions for this case are found to be of the form 

u{t) = A exp (-7r) [Joippr) + PJ^ippr) cos (40)) , (10) 

v{r) = A exp (— 7r) (aJi^ppr) exp + 5 Joippr) exp (3i0)) . (11) 

In analogy to Ref. |T^, we obtain 

E = K-V-^A{a + 5), (12) 



Ea = -Ka-^Ail + ^), (13) 

Ep = Kp-^A{a + 6), (14) 
M = -K<5--^A(1 + A). (15) 

From Eqs. (12) - (15) we get a = 6 = {V - K)/{V2A) and (3 = {1 - V/{K - E))/V2. For 
an impurity with a scattering strength in the unitary hmit the energy of the bound state 
is expected to be very small, ~ 0, and V ^ A. This gives P ~ —1/^/2. On the other 
hand, in the Born limit ~ A and ~ 0, which gives P ~ l/\/2. In Fig. 2(a) and (b), 
we plot |M(r)p and |f(r)p for P = corresponding to weak scattering. In the limit of 

strong impurity scattering, P = — l/y^, the patterns are changed as shown in Fig. 3(a) and 
(b). It is observed that for 2D f-wave superconductors with A(k) = Az cos (20) exp (±Z(/)), 
both |'u(r)p and |f(r)p have a four-fold symmetry. Depending on the impurity scattering 
strength, they extend either in the directions of the Ru-0 bonds or are tilted by an angle of 
45 degrees. The STM imaging of impurity bound states can thus provide a clear signature 
of the order parameter symmetry for the underlying superconductivity. 



III. CONCLUDING REMARKS 



Stimulated by the successful scanning tunneling microscopy imaging of quasi-particle 
bound state wave functions around Zn-impurities in Bi2212 [0,|l^, we have studied the anal- 
ogous patterns of impurity bound states in Sr2Ru04. By applying the appropriate Bogoli- 
ubov - de Gennes equations, the characteristic patterns were distinguished for two proposed 
order parameters: (i) gapped 2D p-wave (or Eu) superconductors with A(k) = Az exp {±i(f)), 
and (ii) gapless 2D f-wave (or Big x Eu) superconductors with A(k) = Az cos (20) exp (±i0). 
While the tunneling conductance patterns of the fully gapped odd-parity 2D p-wave super- 
conductor are featureless along the azimuthal direction, a clear four-fold symmetry in |ii(r)p 
and |f (r)p is predicted for the 2D f-wave superconductor. 

The experimentally observed in-plane anisotropy of the upper critical field Hc2 in 
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Sr2Ru04 is quite small (~ 3%). This may indicate that there is a combination of 2D 
f-wave and p-wave superconductivity in this compound. |[T^ Furthermore, a possible alter- 
native to the plain 2D p-wave superconductivity considered above would be a 3D Aig x Eu 
f-wave order parameter of the form A(k) = Az cos {ckz) exp {±i(f)). The in-plane impurity 
bound state patterns for 2D Eu and 3D Aig x Eu are the same, and a directional probe along 
the fcz-direction would be needed to distinguish between these two cases. 

Our study suggests that the Bogoliubov - de Gennes formalism in the continuum limit 
is very useful in addressing the shape of impurity induced bound states. The corresponding 
bound state wave functions u{r) and v{r) clearly reflect the symmetry of A(r). Therefore 
the imaging of these wave functions provides unique insight into the underlying symmetry 
of the order parameter. A similar study of the vortex state is in progress. 

We thank Kazushige Machida, Manfred Sigrist, and Hyekyung Won for useful discussions. 
K. M. acknowledges the hospitality and the support of the Max-Planck Institut fiir Physik 
komplexer Systeme at Dresden, and S. H. thanks the Zumberge foundation for financial 
support. 



6 



REFERENCES 



[1] S. Nishizaki, Y. Maeno, and Z. Mao, J. Phys. Soc. Jpn. 60, 572 (2000); S. Nishizaki, Y. 
Maeno, and Z. Mao, J. Low Temp. Phys. 117, 1581 (1999). 

[2] K. Ishida, H. Mukuda, Y. Kitaoka, Z. Q. Mao, Y. Mori, and Y. Maeno, Phys. Rev. Lett. 
84, 5387 (2000). 

[3] L Bonalde, B. D. Yanoff, M. B. Salamon, D. J. van Harhngen, E. M. E. Chia, Z. Q. 
Mao, and Y. Maeno, preprint. 

[4] H. Won and K. Maki, Phys. Rev. B 49, 1397 (1994). 

[5] T. M. Rice and M. Sigrist, J. Phys. Cond. Mat. 7, L643 (1995); M. Sigrist, D. Agterberg, 
A. Furusaki, C. Honerkamp, K. K. Ng, T. M. Rice, and M. E. Zhitomirsky, Physica C 
317 - 318, 134 (1999). 

[6] G. M. Luke, Y. Fudamoto, K. M. Kojima, M. L Larkin, J. Merrin, B. Nachumi, Y. J. 
Uemura, Y. Maeno, Z. Q. Mao, Y. Mori, H. Nakamura, and M. Sigrist, Nature 394, 
558 (1998). 

[7] K. Ishida, H. Mukuda, Y. Kitaoka, K. Asayama, Z. Q. Mao, Y. Mori, and Y. Maeno, 
Nature 396, 653 (1998). 

[8] Y. Hasegawa, K. Machida, and M. Ozaki, J. Phys. Soc. Jpn. 69, 336 (2000). 

[9] T. Dahm, H. Won, and K. Maki, |cond-mat/ 0006301 . 

[10] Z. Q. Mao, Y. Maeno, S. Nishizaki, T. Akima, and T. Ishiguro, Phys. Rev. Lett. 24, 
991 (2000). 

[11] S. H. Pan, E. W. Hudson, K. M. Lang, H. Eisaki, S. Uchida, and J. C. Davis, Nature 
403, 746 (2000). 

[12] Y. Morita, M. Kohmoto, and K. Maki, Europhys. Lett. 40, 207 (1997). 



7 



[13] S. Haas and K. Maki, to be published in Phys. Rev. Lett. (|cond-mat /00043"8^ ) . 



[14] For the high-Tc cuprates, it is known that their superconducting coherence length is 
rather short, ^ ~ 10 — 20A. Therefore the impurity induced spatial dependence of the 
gap function A(r) was neglected in Ref. [|I^. In contrast, for Sr2Ru04 ^ ^ 1/im. Hence 
a self-consistent treatment of the spatial variation in A(r) may be needed for a more 
quantitative account of the impurity problem in this material. 

[15] M. Matsumoto and M. Sigrist, J. Phys. Soc. Jpn. 68, 724 (1999); M. Matsumoto and 
M. Sigrist, Physica B 281 - 282, 973 (2000). 

[16] Y. Okuno, M. Matsumoto, and M. Sigrist, J. Phys. Soc. Jpn. 68, 3054 (1999). In this 
reference, impurity induced circular currents were studied. The corresponding wave 
functions are practically equivalent to our solutions. 

[17] I. I. Mazin and D. J. Singh, Phys. Rev. Lett. 79, 733 (1997). 

[18] E. W. Hudson, S. H. Pan, A. K. Gupta, K. W. Ng, and J. C. Davis, Science 285, 88 
(1999). 



8 




-4-2024 -4-2024 



FIG. 1. Spatial variation of the local tunneling conductance, centered at a non- magnetic 
impurity in a 2D p-wave superconductor with A(k) = Az exp (±i^). In the left figure, the dominant 
contribution |n(r)p at the positive bound state resonant frequency Eq is shown. On the right hand 
side, the dominant contribution |f(r)p at the negative bound state resonant frequency —Eq is 
shown. 
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FIG. 2. Spatial variation of the local tunneling conductance, centered at a non-magnetic 
impurity in a 2D f-wave superconductor with A(k) = cos (20) exp {±i<p). On the left hand side, 
the dominant contribution |u(r)p at the positive bound state resonant frequency Eq is shown. 
On the right hand side, the dominant contribution [^(r)^ at the negative bound state resonant 
frequency —Eq is shown. The solution in this figure corresponds to the weak impurity scattering 
limit. 
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FIG. 3. Same as Fig. 2, but for strong impurity scattering. 
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